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Abstract 

We analyze the ability of the three different Liquid Drop Mass (LDM) formulas to describe 
nuclear masses for nuclei in various deformation regions. Separating the 2149 measured 
nuclear species in eight sets with similar quadrupole deformations, we show that the masses 
of prolate deformed nuclei are better described than those of spherical ones. In fact, the 
prolate deformed nuclei are fitted with an RMS smaller than 750 keV, while for spherical 
and semi-magic species the RMS is always larger than 2000 keV. These results are found to 
be independent of pairing. 

The macroscopic sector of the Duflo-Zuker (DZ) mass model reproduces shell effects, 
while most of the deformation dependence is lost and the RMS is larger than in any LDM. 
Adding to the LDM the microscopically motivated DZ master terms introduces the shell 
effects, allowing for a significant reduction in the RMS of the fit but still exhibiting a better 
description of prolate deformed nuclei. The inclusion of shell effects following the Interacting 
Boson Model's ideas produces similar results. 
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I. INTRODUCTION 



When nuclear physicists refer to the description of nuclear masses employing the Liquid 
Drop (LD) model formula, we have a contradictory speech. In fact, sometimes we say: "The 
liquid-drop energy of a spherical nucleus is described by a Bethe-Weizsacker mass formula" 
( 1| citing {2!) or "This is a crude model that does not explain all the properties of the nucleus, 
but does explain the spherical shape of most nuclei" jsj. However, at the same time, we 
say: "The semi-empirical mass formula gives a good approximation for atomic masses and 
several other effects, but does not explain the appearance of magic numbers." jsl, which 
are usually considered the 'more spherical' nuclei, because deformation is associated to the 
quadrupolar interaction between valence protons and neutrons The usual procedure to 
find a phenomenological mass formula which reduce the root mean square (RMS) for the 
actually measured 2149 nuclear species, is to start with the Liquid Drop Mass (LDM) and 
add to it corrections due to deformation and shell-effects. 

The description of nuclear masses in terms of the LDM paved the way to the basic 
understanding of nuclear properties, like the saturation of the nuclear force, the existence 

n 

of pairing and shell effects, and the description of fission and fusion processes [5[. The 
Q-values of different nuclear reactions, obtained from mass differences, must be accurately 
known to allow the description of the astrophysical origin of the elements Q]. Accurate 
theoretical predictions of nuclear masses remain a challenge Q] , sharing the difficulties with 
other quantum many-body calculations, and complicated by the absence of a full theory of 
the nuclear interaction. 

n 

Decades of work have produced microscopic and macroscopic mass formulas [8]. At 
present, the most successful approaches seem to be the microscopic-macroscopic models. 



like the Finite Range Droplet Model (FRDM) [9], its improvements 10|, and the realistic 
Thomas-Fermi (TF) models tU 



12| , the Skyrme and Gog ny Hartee Fock Bogolyubov (HFB) 



and the Duflo-Zuker (DZ) mass formula |l5l-ll7j. They allow for the calculation of 
masses, charge radii, deformations, and in some cases also fission barriers. They all contain 
a macroscopic sector which resembles the LDM formula, and include deformation effects. 
HFB calculations are now able to fit known nuclear masses with deviations competitive with 
the microscopic-macroscopic calculations, while the most precise and robust nuclear mass 
predictions are given by the DZ model [s], [isl, which gives an RMS of 373 keV. 



Efforts for building algebraic nuclear mass formulas inspired in the DZ model success have 
led to detailed analysis of the microscopic building blocks of this model )l9|. which suggested 
new ways to introduce the shell effects through the DZ master terms f20|. Other line of 
thought connected the Interacting Boson Model (IBM) F-spin with the DZ microscopic 



terms 



18 
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22|, which allows for very good fits of the nuclear binding energies when an 



additional one-body Hamiltonian with a large number of parameters is employed [23|. 

Following these works, we perform in this paper an analysis of the interplay between de- 
formations and shell effects for different nuclear mass formulas. We start studying the ability 
of LDM to describe nuclear masses for nuclei in different deformation regions and explicitly 
show that the best fit is obtained for deformed nuclei, extending a previous preliminar work 



24j |. These results are found to be independent of the pairing term, which could be failing 



in the vicinity of closed shell nuclei. 

Analyzing the macroscopic terms of the DZ mass model we show that shell effects are 
reproduced and most of the deformation dependence is lost and the RMS is larger than in 
any LDM. Adding to the LDM the microscopically motivated DZ master terms introduces 
the shell effects, allowing for a significant reduction in the RMS of the fit but still exhibiting 
a better description of prolate deformed nuclei. The inclusion of shell effects through a 
dependence in the number of valence nucleons produces similar results. 

The paper is organized as follows: the fits for three LDM are presented in Sect. 2, and 
the DZ and other microscopic algebraic estimations are discussed in Sect. 3. A meticulous 
analysis of the master terms is performed in Sect. 4, and a comparison with other estimations 
based on IBM is given in Sect. 5. Conclusions are drawn in Sect. 6. 



II. THE FITS FOR THREE LDM FORMULAS 



We have selected three LDM formulas to analyze their ability to fit nuclear masses. 



A. LDMl 



The first one is an improved version of the LDM formula with modified symmetry and 
Coulomb terms, built following a consistent treatment of nuclear bulk and surface effects 



22| . The negative nuclear binding energy is given by 



,2/3 ^ 4T(T + 1) Z{Z-l) A 

where: i) the pairing interaction is given by A = 2, 1, and for even-even, odd-mass and 
odd-odd nuclei, respectively; ii) a correction to the radius of the nucleus is included through 
a modification A in the Coulomb term, A = QZ(i+y-^A^/'-^) ' ™^ symmetry term employs 
4T{T + 1), with T = \N — Z\/2, instead of (A^ — Z)"^ to account for the Wigner energy; and 
iv) the Coulomb interaction is proportional to Z{Z — 1) to avoid the Coulomb interaction 
of a proton with itself. 



B. LDM2 



The second LDM formula is a modified version of the Bethe-Weizsacker one, which in 
corporates exphcitly isospin effects 

6, 



Eldm2 = -ayA + asA^/^ + 



''sym 



I^A + a.^^^{l-Z^'^)-a. 



np 



with isospin asymmetry / = (A^ — Z)/A. The pairing term is taken from 19 | 



(2) 



(3) 



2— I / I :N and Z even, 

I / I :N and Z odd, 

1— I / I :N even, Z odd, and N>Z, 

1— I / I :N odd, Z even, and N<Z, 

1 :N even, Z odd, and N<Z, 

1 :N odd, Z even, and N>Z, 

and the symmetry energy coefficient, including an / correction on the conventional surface- 
symmetry term of LDM to approximately describe the Wigner effect for heavy nuclei, is 
written as 
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C. LDM3 



The third LDM formula is derived from the hquid drop model with the pairing energies 



of the Thomas-Fermi model 



25| 



Eldm3{N, Z) = -a,(l - KP)A + a,(l - ksP)A^'^ + au{l - hP)A^/^ + 



5 roAV3 



(5) 



-fpZ^A - a,,e..Z4/VAV3 + E,,„. 

It includes: i) a first term representing the volume energy corresponding to the saturated 
exchange force and infinite nuclear matter, with PA being the asymmetry energy of the 
Bethe-Weizsacker mass formula; ii) a surface energy corresponding to semi-infinite nuclear 
matter and originated by the deficit of binding energy of the nucleons at the nuclear surface; 
iii) a curvature energy resulting from non-uniform properties which correct the surface energy 
and depends on the mean local curvature; iv) the decrease of binding energy due to the 
Coulomb repulsion; v) a Z'^/A diffuseness correction to the sharp radius Coulomb energy; vi) 
a Z^^'^ /A^^^ charge exchange correction term; vii) the pairing energies Epair of the Thomas- 
Fermi model (taken from Eq. (A. 2) in Ref. 



D. The fits 



The coefficients of the three LD models were selected to minimize the RMS when the pre- 
dicted binding energies BEii^{N, Z) are compared with the experimental ones BE^xpiN, Z), 
reported in AME03 26||, modified so as to include more realistically the electron binding 
energies as explained in Appendix A of Lunney, Pearson and Thibault 8|: 

l2 ^ 1/2 

(6) 



RMS 



E[5Eexp(Ar, Z)-BEty,{N, Z)X 



nucl 



The minimization procedure uses the routine Minuit |2j 

The fits were performed separating the nuclei in nine groups: 



the first one contains all nuclei whose measured masses are reported in AME03 26| . 
which have N, Z > 8, 



the next seven groups contain the nuclei whose quadrupole deformations 62, taken 



form the FRDM 



, lie in the ranges listed in the second and third row of Table HI 



TABLE I: The nine groups of nuclei employed in the present study, their range of quadrupole 
deformation, and their number of nuclei. 



group 


all 


1 2 3 4 5 6 7 


semi-magic 


62 min 


-0.65 


-0.65 -0.11 0.00 0.04 0.12 0.18 0.23 




62 max 


0.65 


-0.11 0.00 0.04 0.12 0.18 0.23 0.65 






2149 


258 252 332 272 307 364 364 


185 



• the last group contains all semi-magic nuclei, having Z = 14, 28, 50, 82 or = 14, 
28, 50, 82 or 126. 

Nnuci is the number of nuclei in each group, listed in the fourth row of Table [II We remark 
that regions 1 to 7 contain approximately the same quantity of nuclei. Notice that group 1 
contains most of the oblate nuclei, that the more spherical nuclei belong to groups 2, 3 and 
4, and that the more prolate deformed nuclei are included in groups 6 and 7. We show in 
Fig. [T]the different regions in the N-Z plane. Semimagic nuclei are displayed along straight 
thick black lines. Around them cluster the nuclei classified as the more spherical, while the 
more prolate-deformed ones form closed regions with many valence protons and neutrons. 

For each LDM equation, nine fits were performed, one for each group of nuclei. In this 
way, nine sets of parameters were obtained, which minimize the RMS of each group of nuclei. 



The values of these parameters can be found in Tables 3, 5 and 7 from Ref. 2J]. We present 
in Table Hi] the RMS obtained for all groups employing these nine sets of parameters. 

TABLE II: RMS (in keV) for each of the nine groups (columns), employing Eqs. ([T]), ([2]) and ([5]), 
respectively. 





all 


1 2 3 4 5 6 7 


semi-magic 


LDMl 


2387 


1313 1676 2063 1746 1053 870 746 


2113 


LDM2 


2374 


1254 1675 2069 1762 1021 838 656 


2056 


LDM3 


2422 


1183 1597 2151 1517 986 819 629 


1967 



These results show that for the three LDM models the masses of prolate deformed nuclei 
can be described with remarkable precision, with an RMS smaller than 750 keV, while the 
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FIG. 1: The seven deformation regions and the semimagic nuclei, shown in the N-Z plane. 

masses of spherical and semi-magic nuclei are those worst described, with RMS larger than 
2000 keV. It completely challenged the authors preconception that the LDM was best suited 
to describe spherical nuclei, and that deformation effects were a crucial element necessary 
to improve the LDM description of nuclear binding energies. As can be read from Table HTl 
deformed regions are very well adjusted by the three LDM formulas, and the regions around 
shell closures are those which represent a challenge. 

Given that the description of pairing effects for nuclei with few valence nucleons would 
require more than a simple parametrization as the one employed here, we have removed odd 
- even effects from our analysis, separating the nuclei in four additional sets, corresponding 
to even A^-even Z (ee), even A^-odd Z (eo), odd A^-even Z (oe) and odd A^-odd Z (oo). The 
number of nuclear species of each type in each set is indicated in Table lllli We repeated the 
fits for each ee, eo, oe and oo set in the nine regions, employing the LDMl mass formula but 
without the pairing term (last term in Eq. ([T])). We arrive to the RMS indicated in Table 



7 



TABLE III: Number of nuclei with even A^-even Z (iV^^^z)' ^ven iV-odd Z (iV^"^/)' odd A^-even Z 
iKlci) and odd A-odd Z (A-,,). 



group 


all 


1 2 3 4 5 6 7 


semi-magic 




576 
534 
536 
503 


62 35 154 53 80 104 88 
67 74 67 79 73 89 85 
69 60 80 62 78 89 98 
60 83 31 78 76 82 93 


95 
40 
50 




IIVI which one more time indicate that prolate deformed nuclei in region 7 are far better 
described than spherical ones, regions 2, 3, 4 and semimagic. It shows that the difficulty 
of the LDM mass formulas to describe spherical nuclei is not at all associated with any 
odd-even effect, or with the way the pairing contribution is parameterized. 

TABLE IV: RMS (in keV) for each of the nine groups (columns) in each ee, eo, oe and oo region, 
employing Eq. ([1]) without pairing. 





all 


1 2 3 4 5 6 7 


semi-magic 


even A-even Z 


2418 


1158 1443 2272 1529 1050 805 786 


2194 


even A-odd Z 


2404 


1387 1834 1698 1726 1036 867 712 


1483 


odd A^-even Z 


2360 


1382 1514 1936 1638 990 880 749 


1871 


odd A^-odd Z 


2404 


1430 1370 1595 1888 1072 913 703 





In what follows we explore the ability of the algebraic extensions of the LDM to include 
shell effects in the description of nuclear masses. 

III. THE DUFLO-ZUKER INSPIRED MASS MODELS 

In this section two different mass models, both based on Duflo-Zuker ideas, are discussed 
in detail. The master terms are introduced in two alternative way in order to describe the 
shell effects in mass models, and their correlation with the nuclear deformation is analyzed. 
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A. DZl 



The macroscopic sector of the simplest version of the DZ mass model, contains six terms 



leading asymptotically to a LD form 



19| (see section C below) 



E 



DZl 



ai{M"^ + S)-a2 



M 



HO 



a^Vc - a^Vx + 05^^.5 + ae^p. 



(7) 



The master M^'^ term is 



M 



HO 



1 

2p 



E 



DHO 

V 



or 



(8) 



where the sums run over all occupied proton and neutron orbitals up to the Fermi level, 



D^o ^ + + 2) is the de generacy of the major harmonic-oscillator (HO) shell of 
principal quantum number p, and n^{p), tLt^^p) are the number of neutrons and protons, 
respectively, in the HO shell p. 

It is relevant to mention that this form describes the dominant contribution of the 
monopole part of the nuclear Hamiltonian. In order to change the HO closures {aX N, Z = 8, 
20, 40, 70, ■ ■ ■) into the observed extruder-intruder (EI) ones at N, Z = 14, 28, 50, 82 and 
126, the 5* operator proposed by Duflo, given in Eq. (18) of Ref. 19|], is employed. 

The scaling factor is 



a) 



212 



p = A'/' 

The Coulomb, asymmetry and surface asymmetry terms are [19(], respectively: 
-Z(Z- 1) + 0.76[Z(Z- 1)2/3] r 



(9) 



Vc 
Vt 
Vts 



^1/3 



4T(T + 1) 

4T(T+1) 4T(r-|) 



The pairing term Vp is the same employed in LDM2. 



(10) 



B. DZ2 



As there exists some uncertainty in the parameterization of the monopole part of the 
nuclear Hamiltonian 19|, in Ref. 20| a modified master term M was proposed, which 



9 



directly includes the EI shell closures. It is built using an expression similar to ([8]), but with 
the index p now referring to the EI major shell with degeneracies Dp — )■ D^'^ + 2. 

MA , , 

-Edz2 = axMA - a2 asVc - a^VT + a^Vrs + a^Vp, (11) 

P 

with the new master term: 

MA = -{el + el), e^. = E e,. = E (12) 

where Dp^^ = (^^,7^ + l)(Pi^,7r + 2) +2 contains the HO — )■ EI information and ny{p^), Ut^^Pt,) 
are the number of neutrons and protons in their respective EI shell. 



C. The master term 

The master term MA contains the microscopic information about shell closures. We 
show in Fig. [2] its behavior as a function of A. It scales as a linear function in A. Therefore, 
we can interpreter this term as a volume one and, consequently, the term MA/p can be 
associated to a surface one. This allows us to consider the macroscopic sector of the DZ 
mass formula as an extension of the LDM. 

The building blocks of the master term MA in Eq. f[T^ are eii, and Ci^. In Fig. [3] the 
behavior of function of is presented. The thin black line represents the asymptotic 



form eiiy ~ ^^^^ ^ 20 1. The plot confirms that Ci,^ scales as A^^^^, with minor differences for 



small A^. To visualize this deviation, the difference between cii, and this asymptotic behavior 
is presented in Fig. H] (black triangles). The microscopic elements introduced through ei^ are 
made evident in this plot. When the asymptotic behavior is removed, shell effects emerge, 
with well defined peaks at shell closures. However, there is a remnant continuous increase 
with N which has not a linear dependence. The description of the asymptotic form of ei^, 
can be improved by including terms with different powers of N^^^. The differences between 
Cii, and these asymptotic forms are displayed with different symbols in Fig. HI From these 
results we conclude that a better description can be reached if we adopt the new asymptotic 
form 

ei^,asym = "0.90892 + 0.54259A^^/^ + 0.98851Ar2/^ ^ 0.0018Ar. (13) 
A similar analysis was performed around Fig. 3 of Ref . jlj] . 



in 




FIG. 2: Master term from DZ model as function of A: gray squares indicate the results calculated 
with Eq. ()12p and solid black line is a linear fit of them. 

In the construction of the master term MA, Eq. ( fT2|) . we employ e^^,. It is compared 
with its asymptotic form ej^^^^^ in Fig. O The advantages of using the Eq. ( 1T5]) instead 
of the simpler ^^^^ ^ expression can be clearly observed: we obtain a global leveling besides 
the good description of closed shell effects. 

Thus the master term MA has the asymptotic behavior described by 



- (e^ + (U) 

y lu,asym ' Ivr, asymy ' \ / 



D. The fits 



Employing the two versions of the macroscopic sector of the DZ mass model, Eqs. ([7]) 
and ffTTj) . we have fitted the nuclear masses in the nine regions described above. 

The fitting procedure in the different deformation regions leads to the RMS displayed in 
Table |Vl These results show that the ability of both models to describe masses of nuclei in 
spherical, prolate and semi-magic groups are now comparable. The global RMS are larger 
than those obtained with the LDM formulas. It is hard to find any correlation between the 

11 
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FIG. 3: Dependence of ei^ with N. 
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FIG. 4: Difference between eii, and different approximations to describe its asymptotic behavior 
as a function of N. 

RMS and the regions with different deformations. This fact is also reflected in the values 



19 




FIG. 5: Difference between ef^, and its asymptotic behavior for the approximation eii, ~ -^^ — ^ — 
(gray hne) and for Eq. (jl3p (black line). 

TABLE V: RMS (in keV) for each of the nine groups, using Eqs. ^ and ([TT]) . 





all 


1 2 3 4 5 6 7 


semi-magic 


DZl 
DZ2 


2852 
3443 


994 1969 1557 1392 2237 2562 1529 
1425 1544 2167 1717 1729 2047 2107 


1392 
1973 



of the fitted parameters, presented in Table IVIII and IVIIII of the Appendix, because some of 
them (mainly and ag) vary noticeably from one region to the other. 

IV. ADDING SHELL CORRECTIONS TO THE LDM 
A. LDMl plus DZ 

It is possible to combine in a single mass formula both the ability of the LDM to describe 
the masses of deformed nuclei, and of DZ to include shell effects. To this end, we will 
work with a LDM, adding the shell effects present in the DZ master term, but with their 
asymptotic behavior removed, because it is basically contained in the volume and surface 

1.'^ 



terms of the LDM. We have selected as the starting point the LDMl, which has a global 
RMS smaller than those of DZl and DZ2. 

Based in our previous analysis, the shell effects will be introduced by a "volume" and 
"surface" shell corrections defined as MA — MAasym and {MA — MAasym)/p, respectively. 
They are plotted in Fig. [6] function of A. 
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FIG. 6: Volume {MA — MAasym) and surface {{MA — MAasym) / p) terms as function of A. 

They resemble the shell effects not included in the LDM. To help the visualization of the 
shell effects introduced by these differences, we show in Fig. [7] the same terms as function 
of A^, Z and A only for semi-magic nuclei, where we can observe one more time the clear 
peaks at the magic numbers. 

As explained above, we have added the shell effects to the LDMl to construct a new mass 
formula: 



ElDM+DZ = LDMl + ayol{MA - MAasym) + asurf{MA - MAasym) /p- 



(15) 



B. Comparison with other estimations of shell effects 



Inspired in the F-spin symmetry of the IBM's 



21 



22|, shell effects have been recently 



introduced in extensions of the LDM by adding two terms, linear and quadratic in the 
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FIG. 7: Volume (MA — MAasym) and surface {{MA — MAasym)/p) terms for semi-magic nuclei. 



number of valence protons, and neutrons, 
the nuclear binding energies is 



21 



, |22| . The proposed expression for 



EiDM+vai = LDMl + 6i(n„ + z^) + h2{ny + z^f. 



(16) 



C. The fits 



Repeating the analysis for the same nine groups, we obtain the results shown in Tables 
IVII for the RMS employing the two models. 

TABLE VI: RMS (in keV) for each of the nine groups, using Eqs. (jl5p and (jl6p . 





all 


1 2 3 4 5 6 7 


semi-magic 


LDM+DZ 


1407 


668 907 1026 755 784 791 647 


952 


LDM+val 


1075 


796 981 1006 828 711 836 615 


1037 



As expected, the improvements in the inclusion of shell effects reduce the global RMS from 
2387 keV to 1407 (1075) keV when the LDM+DZ (LDM+val) model is used. The results 
still show a visible tendency to describe better the deformed than spherical nuclei. Besides, 
the results obtained with both formulas look very similar, with a smaller global RMS in the 
valence model and some advantage of the LDM+DZ model to describe semi-magic nuclei. 

The parameters set obtained in both fits are presented in the Appendix, in Tables IIXI 
and 1X1 The relative stability of the parameters hi and 62 in the LDM+val model, shown 

in 



in the last two columns of Table IXl is behind the comparatively small global RMS. On the 
other hand, the shell surface and volume coefficients ag^rf and a^oi of the LDM+DZ model 
listed in the last two columns of Table IIXI vary both in magnitude and in sign from one 
deformation region to another. Being a limitation for a good global fit, it offers at the same 
time the opportunity to relate these parameters with the deformation, a challenge which is 
left for future work. 

Finally, to complete the comparison, we exhibit in Fig. |8]the difference between exper- 
imental data (Bexp) and theoretical results (Bth) for both models. The two models which 
include shell effects display a clear reduction in the differences as compared with the Liquid 
Drop Model LDMl. As decades of work in nuclear masses have thought us, including these 
residual effects is a fairly non-trivial task. 




FIG. 8: Comparison between LDMl, LDM+DZ and LDM+val fits. 



V. CONCLUSIONS 

We have explored the ability of several Liquid Drop Models to describe nuclear binding 
energies in different deformation regions, and we have performed a similar analysis employing 
more elaborated models which include shell effects through microscopic terms. We selected 
for these microscopic corrections the Duflo-Zuker type models and those inspired in the 
Interacting Boson Model F-spin. 

We have shown that the LDM is best suited to describe the masses of prolate deformed 
nuclei than of spherical ones, while these deformation effects are washed out employing the 

Ifi 



macroscopic sector of the DZ mass formulas. 

We have paid special attention to the DZ master term, which can be used to construct 
volume and surface terms which describe shell effects. Adding these microscopic information 
to the LDM formula, we obtain an eight parameter fit of nuclear masses, with a global 
RMS of 1407 keV. The strong dependence of the parameters of the new shell volume and 
surface terms of the deformation regions could open the possibility to relate shell effects and 
deformation in an more elaborated way. 

Acknowledgments 

C.B. and A.M. are fellows of the CONICET, CCT La Plata (Argentina). JGH thanks A. 
Aprahamian, O. Civitarese, S. Pittel and P. van Isacker for their valuable comments. This 
work was supported in part by the Agenda Nacional de Promocion Cientifica y Tecnologica 
through PICT-2007-00861, Conacyt, Mexico, and DGAPA, UNAM. 



[1] N. Wang, M. Liu and X. Wu, Phys. Rev. C81 (2010) 044322 . 
[2] H.A. Bethe and R.F. Bacher, Rev. Mod. Physics 8 (1936) 82-229. 



[3] http://en.wikipedia.org/wiki/Semi-empirical mass formula. 



[4] P. Federman and S. Pittel, Phys. Rev. C20, (1979) 820. 

[5] A. Bohr and B.R. Mottelson, Nuclear Structure v. 1, (World Scientific, Singapore, 1998). 
[6] C.E. Rolfs and W.S. Rodney, Cauldrons in the Cosmos (University of Chicago Press, Chicago, 
1988). 

[7] Klaus Blaum, Phys. Rep. 425 (2006) 1. 

[8] D. Lunney J.M. Pearson, C. Thibault, Rev. Mod. Phys. 75 (2003) 1021. 

[9] P. Moller, J.R. Nix, W.D. Myers, W.J. Swiatecki, At. Data Nucl. Data Tables 59 (1995) 185. 
[10] N. Wang, Z. Liang, M. Liu and X. Wu, Phys. Rev. C 82 (2010) 044304. 
[11] W.D. Myers and W.J. Swiatecki, Nucl. Phys. A 601 (1996) 141. 
[12] K. Pomorski and J. Dudek, Phys. Rev. C67, (2003) 044316. 

[13] S. Goriely, F. Tondeur, J.M. Pearson, Atom. Data Nucl. Data Tables 77 (2001) 311; S. Goriely 
M. Samyn, J.M. Pearson, Phys. Rev. C75 (2007) 064312. 



17 



[14] S. Goriely, N. Chamel, J.M. Pearson, Phys. Rev. Lett. 102 (2009) 152503; S. Goriely, S. 

Hilaire, M. Girod, adn S. Peru, Phys. Rev. Lett. 102 (2009) 242501. 
[15] J. Duflo, Nucl. Phys. A576 (1994) 29. 
[16] A.P. Zuker, Nucl. Phys. A576 (1994) 65. 
[17] J. Duflo and A.P. Zuker, Phys. Rev. C52 R23 (1995) R23. 

[18] J. Mendoza-Temis, I. Morales, J. Barea, A. Frank, J.G. Hirsch, J.C. Lopez- Vieyra, P. Van 

Isacker and V. Velazquez, Nucl. Phys. A812 (2008) 28. 
[19] J. Mendoza-Temis, J.G. Hirsch and A.P. Zuker, Nucl. Phys. A843 (2010) 14. 
[20] J.G. Hirsch and J. Mendoza-Temis, J. Phys. G: Nucl. Part. Phys. 37 (2010) 064029. 
[21] A.E.L. Dieperink, P. Van Isacker, Eur. Phys. J. A 32 (2007) 11. 
[22] A.E.L. Dieperink, P. Van Isacker, Eur. Phys. J. A 42 (2009) 269. 
[23] G. Gangopadhyay, Int. J. Mod. Phys. E 20 (2011) 179. 

[24] J.G. Hirsch, C. Barbero and A. Mariano, larXiv: 1108. 0707] J. Phys.: Conference Series (in 
press) . 

[25] G. Royer, M. Guilbaud and A. Onillon, Nucl. Phys. A847 (2010) 24. 
[26] G. Audi, A.H. Wapstra, C. Thibault, Nucl. Phys. A 729 (2003) 337. 

[27] F. James, Minuit: Function Minimization and Error Analysis Reference Manual, Version 94.1, 
CERN (1994); |http : / / wwwasdoc . web . cer n . ch / wwwasdoc / minuit /minmain . html 

[28] P. Moller, J. R. Nix, W. D. Myers, and W. J. Swiatecki, Atomic Data Nucl. Data Tables 59 
(1995) 185. 



18 



Appendix: Parameters of the fits 



In this appendix we present the parameters found to provide the best fits of the nuclear 
masses in the nine regions for the mass equations DZl (Eq. ([7])), DZ2 (Eq. ( ITTl) ). LDMl+DZ 
(Eq. (US])) and LDMl+val (Eq. (^). 



TABLE VII: Sets of parameters (in keV) which minimize the RMS for the nine groups of nuclei, 
employing DZl, Eq. ([7|). 



seti 
set2 
sets 
set4 
sets 
sete 
set 7 

set cfim^ 



ai 



0.2 03 



04 



05 «6 



17.531 
17.458 
17.760 
17.788 
17.975 
17.894 
17.172 
17.379 
17.993 



15.433 0.694 
15.204 0.691 
16.283 0.706 
16.329 0.708 
16.971 0.716 
16.593 0.715 
14.375 0.667 
15.037 0.682 
17.054 0.719 



36.157 
33.959 
38.212 
39.263 
39.733 
38.553 
35.674 
33.507 
38.218 



48.425 5.154 
39.860 5.857 
57.641 6.167 
61.183 5.807 
61.192 6.056 
56.304 5.997 
47.421 4.684 
38.541 5.206 
54.671 7.428 



IP 



TABLE VIII: Sets of parameters (in keV) which minimize the RMS for the nine groups of nuclei, 
employing DZ2, Eq. (fTTI) . 



setaii 
seti 
set 2 
sets 
set4 
sets 
setg 
set 7 

set ef=m7 



ai 
15.707 
15.798 
15.801 
15.535 
16.248 
16.116 
15.685 
15.381 
16.204 



02 as 
16.155 0.611 
16.383 0.619 
16.357 0.620 
15.684 0.599 
17.924 0.646 
17.509 0.635 
16.315 0.598 
15.076 0.589 
18.026 0.638 



04 
29.800 
29.892 
30.997 
31.938 
31.604 
32.334 
31.551 
25.663 
32.570 



05 06 

34.307 7.337 

33.626 6.882 

37.886 3.960 

45.543 6.855 

36.711 8.869 

42.394 6.581 

42.802 6.118 

18.329 7.087 

43.255 11.128 



TABLE IX: Sets of parameters (in keV) which minimize the RMS for the nine groups of nuclei, 
using LDMl+DZ, Eq. 1^. 



Os Oc CLp y 



setaii 
seti 
set2 
seta 
set4 
sets 
set6 
sety 
set .00771,1 



15.817 18.422 0.705 5.917 29.740 2.298 2.700 12.526 
15.727 18.142 0.700 5.893 28.721 2.090 17.792 -47.797 
15.889 18.646 0.710 4.980 32.120 2.937 18.350 -48.154 
15.805 18.409 0.703 5.712 33.574 3.352 20.536 -52.987 
16.026 19.043 0.719 5.624 32.050 2.577 26.133 -80.211 
15.911 18.717 0.711 5.731 31.079 2.532 1.962 13.575 
15.703 18.103 0.696 5.297 29.991 2.368 -13.570 67.354 
15.515 17.641 0.682 5.267 27.509 1.894 -10.468 47.078 
15.992 19.008 0.717 7.551 31.819 2.616 22.776 -59.556 



9.0 



TABLE X: Sets of parameters (in keV) which minimize the RMS for the nine groups of nuclei 
using LDMl+val, Eq. ([TH]) . 



an Or. 



y 



15.887 
15.805 
15.870 
15.793 
16.075 
15.929 
15.872 
15.678 
16.091 



18.232 
18.067 
18.038 
17.927 
18.789 
18.397 
18.456 
17.850 
18.824 



0.712 5.651 
0.707 5.882 
0.715 4.555 
0.704 5.115 
0.725 5.086 
0.715 5.715 
0.708 5.795 
0.695 5.083 
0.726 6.796 



31.242 
29.052 
29.963 
32.858 
31.544 
31.276 
30.722 
29.023 
31.392 



2.565 0.771 
2.073 0.606 
2.224 0.848 
3.138 0.780 
2.401 0.730 
2.529 0.781 
2.435 0.380 
2.146 0.456 
2.418 0. 



-0.013 
-0.008 
-0.017 
-0.014 
-0.012 
-0.017 
-0.006 
-0.009 
-0.017 
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